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Abstract 

The piano soundboard transforms the string vibration into sound and therefore, its 
vibrations are of primary importance for the sound characteristics of the instru- 
ment. An original vibro-acoustical method is presented to isolate the soundboard 
nonlinearity from that of the exciting device (here: a loudspeaker) and to measure 
it. The nonlinear part of the soundboard response to an external excitation is quan- 
titatively estimated for the first time, at ~ — 40 dB below the linear part at the jf 
nuance. Given this essentially linear response, a modal identification is performed 
up to 3 kHz by means of a novel high resolution modal analysis technique (Ege et 
al, High-resolution modal analysis, JSV, 325(4-5), 2009). Modal dampings (which, 
so far, were unknown for the piano in this frequency range) are determined in the 
mid-frequency domain where FFT-based methods fail to evaluate them with an 
acceptable precision. They turn out to be close to those imposed by wood. A finite- 
element modelling of the soundboard is also presented. The low-order modal shapes 
and the comparison between the corresponding experimental and numerical modal 
frequencies suggest that the boundary conditions can be considered as blocked, ex- 
cept at very low frequencies. The frequency-dependency of the estimated modal 
densities and the observation of modal shapes reveal two well-separated regimes. 
Below ~ 1 kHz, the soundboard vibrates more or less like a homogeneous plate. 
Above that limit, the structural waves are confined by ribs, as already noticed by 
several authors, and localised in restricted areas (one or a few inter-rib spaces), 
presumably due to a slightly irregular spacing of the ribs across the soundboard. 
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1 Introduction 

Since tlie strings of a piano are too tliin to radiate sound, tliis function is en- 
sured by the soundboard, a complex plate made of many wooden parts. The 
sound of the instrument is therefore largely determined by the vibrational 
characteristics of the soundboard. As shown in Figs. [T] and [7], the main ele- 
ment is a large and thin panel, made out of glued strips (usually in spruce). 
We define the x-direction as the grain direction of this panel's wood. In the y- 
direction are glued a series of parallel, nearly equidistant stiffeners, called ribs 
and also made out of spruce (sometimes sugar pine). On the opposite face are 
glued the bridges: one short and one long thick bars running approximately in 
the x-direction, slightly curved, made out of maple, on which the strings are 
attached. The overall shape of this ensemble depends on the type of piano: 
nearly rectangular for upright pianos, the shape of a half round-and-high hat 
for grand pianos. The width of the soundboard is more or less 140 cm, corre- 
sponding to that of the keyboard. The height or length ranges from more or 
less 60 cm for small uprights to more than 2 m for some concert grands. The 
panel thickness h is between 6 and 10 mm, the inter-rib distance p ranges from 
10 to 18 cm (depending on pianos). The treble bridge is usually 3 to 4 cm high 
and 2.5 to 3.5 cm wide. The bass bridge is about 2-3 cm higher than the treble 
one so that the bass strings can be strung over the lower mid-range ones. The 
dimensions of the cross section of the ribs are about 25 mm in width with 
height varying from 15 to 25 mm (thinner toward the treble region). 

The functioning of the piano goes schematically as follows: once a key has 
been struck, a hammer strikes one, two, or three tuned strings and goes back 
to its rest position, leaving the strings vibrating freely. The bridge represents 
a nearly fixed end for the strings so that energy stored initially in the strings is 
slowly transferred from the strings to the soundboard (a significant quantity 
of energy being also dissipated inside the strings themselves), partly dissi- 
pated inside wood and partly radiated acoustically. Since the decay time is 
several orders of magnitude larger than the periods of vibration, it can be 
considered that the soundboard is put into a forced motion by the strings, 
at frequencies that have no relationship with the resonance frequencies of the 
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Figure 1. Soundboard of a grand-piano (http://www.lindebladpiano.com). Left: up- 
per face, with bridges (bars where strings are attached) visible. Right: lower face, 
with ribs. Note that ribs are closer one to each other in the treble range of the 
instrument (lower part of the picture). 

soundboard itself. Since the string and the soundboard can be considered as 
almost dynamically uncoupled, it makes sense to analyse the dynamics of the 
soundboard in terms of its normal modes. 



The present article is focused on some main features of the vibration regimes 
of the soundboard of an upright piano, as observed in playing condition: lin- 
earity, modal density and modal damping. The literature on the vibrations of 
the soundboard has been recently reviewed in [1]. Curiously, the literature is 
almost mute on the first point, even though linearity is a requisite for the usual 
representations of the dynamics of the piano soundboard - modes, mechani- 
cal impedance or mobility - to be physically meaningful quantities. Section [2] 
presents the first, to the best of our knowledge, quantitative evaluation of 
the (non) linearity of the soundboard vibration. Modal analyses are presented 
next: experimental in Sec. [3] and numerical, by means of a finite-element mod- 
elling of the soundboard in Sec. |4j The frequency range under scope in this 
article [0 - 3] kHz is much wider than in most previous experimental studies 
devoted to modal parameters of the piano. This could be achieved thanks to 
the use of a recently published high-resolution modal analysis technique [2j. 
Compared to techniques based on the Fourier transform, it avoids the cus- 
tomary compromise in time-frequency resolution and thus, gives access to an 
extended frequency-range. Since the timbre of piano depends highly on the 
relative decay-times of the components of each note, energy dissipation repre- 
sents an important dynamical parameter which is accessed here via the modal 
dampings of the soundboard. The evaluation of modal dampings could be per- 



formed here up to several kHz for the first time in piano studies (Sec. 5.1). 
The results are presented and discussed in Sec. [5] in terms of modal shapes in 
the low-frequency range and in terms of modal densities and dampings up to 
3 kHz. Analyses are done in the spectral domain throughout the whole paper 
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and physical quantities are thus considered as imphcit or exphcit functions of 
the frequency /. 



2 A linear behaviour ? 

As recalled above, linearity of the piano must be guaranteed if the usual rep- 
resentations (modes, impedance) are to be used. Moreover, people have ex- 
pressed the opinion that nonlinearities could be significant in the sound of the 
instrument. Both questions reduce to a quantitative aspect (how small?) but 
with the same reference: a 1% distortion rate (-40 dB) may be considered as 
sufficient small for using linear physical concepts but would be considered as 
mediocre for audio equipment. 

2.1 Past studies 

Nonlinear phenomena such as jump phenomenon, hysteresis or internal res- 
onance appear when the transverse vibration of a bi-dimensional structure 
exceeds amplitudes in the order of magnitude of its thickness j3]. In the case 
of the piano, the soundboard transverse motion measured at the bridge re- 
mains in a smaller range, even when the piano is played ff in the lower side 
of the keyboard. Askenfelt and Jansson [4J report maximum values of the dis- 
placement at the bridge of 6 • 10^® m in the frequency range [80-300] Hz 
(Fig. |2|. This maximum value is less than 10~^ times the board thickness. We 
can therefore expect that, to a high level of approximation, the vibration of 
the soundboard is linear. 

Very few experimental studies have been carried out concerning the linearity 
of the piano soundboard vibrations. The most convincing work seems to be 
that by Hundley et al. [5|. Their motivation was in fact to eliminate nonlin- 
earity as one possible cause of the multiple time-decay of piano tones. There- 
fore, these authors studied the nonlinearity in the so-called string-to-bridge- 
to-soundboard-to-air path. A direct measure of linearity was obtained by using 
a magnetic driver close to a string triplet (note B3) and driving the string(s) 
at a fundamental frequency of 250 Hz. The curve shown in Fig. |3](a) repre- 
sents the sound pressure level near the piano in an anechoic environment for 
different input levels of the driver. Proportionality between the sound pressure 
level and the excitation level is excellent up to a SPL of 90 dB, with the slope 
of the curve being very close to one 0.995). However, our opinion is that no 
quantitative evaluation of the nonlinearity can be withdrawn from this curve. 
Since the main response is linear, nonlinearities manifest themselves only at 
high vibration levels so that the major part of this curve do not carry useful 
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Figure 2. Vibration levels at the bridge of a grand piano when played pp (dash-dotted 
line), mf (dotted line) and ff (solid line with • marks) for the notes C2 to B5 
(fundamental frequencies ~ 60 to 950 Hz), according to Askenfelt and Jansson [4j. 
The region below the threshold of vibration sensation at the fingertips (contact area 
28 mm^) is shown shaded. 

information. In order to extract an order-of-magnitude for the nonlinear part, 
the experimental precision at the highest point must be at least one order of 
magnitude higher: 60 dB for measuring a 1% distortion, for example. This is 
clearly not the case in this measurement. 



In a second step, the authors compared the tone level recordings for which 
the key was actuated by different weights (from 100 grams to 800 grams). 
This experiment was repeated for a large number of keys and no evidence of 
dependence between the decay rate and the blow force was found: the curves 
were almost identical (see Fig|3](b) for an example of this measurement on 
note C5). Contrary to the first one, this second experiment is qualitative in 
nature but yields here a null result. 



2.2 Method to estimate the piano soundboard nonlinearities 



We aim at estimating experimentally the order of magnitude for the ratio 
between the nonlinear and the linear parts of the soundboard response - or 
distortion rate - when it is excited by a given string force on the bridge. 
Since replacing the string excitation by a mechanical shaker raises all sorts 
of experimental problems (attachment, feedback of the piano response on the 
shaker excitation, control of the side forces, etc.), we have preferred to excite 
the piano acoustically (with a very powerful electrodynamic loudspeaker) and 
to evaluate its vibratory response, as shown in Fig. |4j 

We assume that, for a given vibratory level of the soundboard, the distortion 
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Figure 3. Linearity of response for the string-to-bridge-to-soundboard-to-air path, 
according to Hundley et al. [5]. (a) with an electromagnetically excited string (after 
Fig. 5). (b) with different weights actuating the C5 key (after Fig. 6). 




Figure 4. Acoustical excitation of the piano placed in a pseudo-anechoic room. The 
acceleration of the board is measured at points Ai, A2 and A5 shown in Fig. [7] 

ratio is the same, at least in order of magnitude, whether the piano is excited 
acoustically or mechanically. If the piano soundboard were a point, the for- 
mer and the latter situations would represent reciprocal experiments and the 
assumption would be exactly true. 

Whether the electrodynamic exciter is acoustical (loudspeaker) or mechani- 
cal (shaker), its own nonlinear contribution to the overall response cannot in 
general be neglected, compared to that of the piano. If the driver response is 
independent from the piano response (which is true to a large extent in the 
acoustical case, for a sufficiently powerful loudspeaker), the (free-of- feedback) 
situation can be represented as in Fig. |5j where the loudspeaker is represented 
by the system Sp, the piano soundboard by Sq and their association by Sh- 
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Figure 5. A chain of two nonlinear systems models the piano soundboard Sq excited 
by the acoustical field created by an electrodynamic loudspeaker Sp- The electrical 
driving signal is X{f), Y{f) is the acoustical signal (see text for discussion) and 
Z{f) is the piano vibratory response. 

Experimentally, only the input electrical signal of the loudspeaker X{f) can 
be considered as fully controlled. Since both the exciter and the soundboard 
are expected to be slightly nonlinear, a direct characterisation of the piano 
soundboard appears to be very difficult. The method which is presented here 
derives the distortion rate of Sq from measurements performed on Sp and Sh- 

The signal Z{f) denotes the acceleration of the soundboard, possibly at various 
different locations. For the sake of nonlinearity estimation, it is assumed that 
the acoustical field created by the loudspeaker can be represented by a scalar 
value Y{f). Moreover, it is assumed that Y{f) can be estimated by removing 
the piano and setting a microphone at the place of the soundboard. Again, the 
numerous corresponding approximations (change of the acoustical field with or 
without the piano, difference between the acoustical field and its measurement 
in one point, etc.) are supposed to be correct only as far as the order of 
magnitude of the nonlinearities is concerned. In other words, we consider that 
the piano, equipped with accelerometers, behaves, as far as nonlinearities are 
concerned, like a slightly nonlinear (and localised) microphone. It is explained 
below how Sp and Sh have been characterised and how one can derive the 
distortion rate of So (the piano soundboard). 

The outputs of the systems Sp, Sg and Sh can generally be decomposed in 
their linear and nonlinear parts as follows: 

Sp: Y{f)^F{f)X{f)[l + CsAf)] (1) 
Sg: Zif) = G{f)Y{f)[l + CsaU-)] (2) 
Sh: Z{f)^H{f)X{f)[l + CsAf)] (3) 

where Csp{f), Cs^{f) and Csffif) are the distortion rates of Sp, Sg and Sh 
respectively. 

It is assumed that the nonlinearities are mathematically weak: the relation- 
ship between their output s{t) and their input e{t) can be expressed ana- 
lytically by Volterra series [6|7] (see Appendix [A]). Let the Volterra kernels 
{Fkifi,-,fk)}km*, {Gfc(/i,...,/fc)}fceN* and {Hkifi, fk)}keN' describe the 
systems Sp, Sg et, Sh- It is shown in Appendix |A] that, in the chain-case 
presented in Fig. [5} the Volterra kernels of Sh are given analytically by func- 
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tions of the Volterra kernels of Sp and Sq- For k = 1, one obtains: 



F^i^) = Fif,) Gi(/i) = G(/i) H,{h) = H{U) (4) 
H{f) = F{f)G{f) (5) 



This result can be extended to a chain of more than two systems and proves 
the intuitive result that the linear part of the response of a chain of weakly 
nonlinear systems is the product of the linear parts (transfer functions) of each 
system composing the chain. 

One assumes further that Sp and Sh can be modelled as cascades of Ham- 
merstein models. Cascade of Hammerstein models constitute an interesting 
subclass of Volterra systems for which parameters can be estimated experi- 
mentally pi9] . as explained in Appendix [B] The purpose here is to derive an 
estimation of Sq (the piano soundboard excited by the acoustical field) based 
on experimental estimations oi Sp (the loudspeaker excited by the driving 
electrical signal) and Sh (the piano combined with the loudspeaker, excited 
by the electrical driving signal). 

Using Eqs. ([s]), ^ and multiplying by Fi(f), the distortion rate Csa{f) of 
the piano is given by: 

F{f)Z{f)-HU-)Y{f) 
Csaif)- j^^j^^ (6) 

where all quantities in the second member of Eq. ^ can be measured as 
described in the next section. This method solves the problem of isolating the 
nonlinearities of a system from those of its exciter. 



2.3 Experimental implementation and results 



An upright piano (Atlas brand) with a rectangular soundboard (dimensions: 
0.91 m X 1.39 m X 8 mm) was used for the experiments. Since the geometries 
of piano soundboards are rather similar (thickness of the wood panels, height 
of the rib, width of the soundboard, width/length of the soundboard varying 
in a 1 to 2 ratio between grands and uprights), the result obtained here can be 
expected to be similar, in order of magnitude, on other pianos. The piano was 
put in a pseudo-anechoic room (anechoic walls and ceiling, ordinary ground). 
It was tuned normally, with strings muted by strips of foam (or woven in two or 
three places) between them. The electrical excitation X{f ) of the loudspeaker 
(Bose - 802 Series II) was an exponential swept-sine [50-4000] Hz (40 kHz 
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sampling frequency, T = 26 s duration). The amplitude of the loudspeaker 
was adjusted so that the displacement of the soundboard did correspond to 
the ff playing level: ^ 10^^ m at 370 Hz in this case. 

The acoustic response Y{f) was measured by a microphone (pre-polarised 
pressure-field 1/2" - Briiel & Kj^r 4947) taking place of the piano. It exhibits 
some distortion which may safely be attributed to the loudspeaker rather 
than to the microphone (the typical distortion rate of the microphone is 3% 
at 160 dB SPL, it becomes totally negligible at a SPL less than 100 dB, 
compared to that of a loudspeaker, around 2-3% at 100 dB). 

The motion Z{f ) of the soundboard was measured with three accelerometers 
(Briiel & Kjasr 4393) put at the locations marked by Ai, A2 and A5 in Fig. [7j 

The distortion rate Cs^if) appears to be generally slightly larger than Csp{f)- 
This is consistent with presumably small soundboard nonlinearities and with 
the expectation that nonlinearities do not compensate each other: nonlinearity 
is expected to increase along the chain. It was also observed that the second 
order distortion is significantly larger than the third-order distortion and much 
larger than the other orders, for both Sp and Sh- In other words, the second- 
order (and principal) nonlinearity of the loudspeaker seems to extend to the 
overall "loudspeaker-soundboard" system. Since the soundboard is nearly a 
flat structure, its intrinsic nonlinearity is expected to be geometric, and of third 
order. However, the present method does not allow to identify the different 
nonlinearity orders in Csaif)- 

In some restricted frequency ranges and for points A2 and A5, it was observed 
that Csp{f) was slightly larger than Cs^lf). Presumably, occurrences of this 
anomalous situation can be attributed to the failure of the hypothesis that 
the piano behaves like a localised system. 

The distortion rates of the piano soundboard Csdf) at points Ai, A2 and A5 
are shown in Fig. [6j The missing parts correspond to the anomalous situation 
described above. The apparent increase in nonlinearity near 4 kHz is probably 
an artefact of the method since the quality of the reconstruction of the nonlin- 
ear impulse responses is bad near the lower and upper bounds of the explored 
frequency range (50-4000 Hz in the present case). The three curves differ by 
not more than an order of magnitude, which is consistent with the fact that 
the level of the soundboard motion differs, but not widely, from place to place 
and therefore, that the responses in different locations do not carry the same 
amount of nonlinearity. Conversely, the fact that the three curves are roughly 
similar justifies, at least for the purpose of finding an order-of-magnitude, the 
approximations made above. 

Altogether, the nonlinear part of the piano response appears to be contained 
within -30 to -50 dB. The order of magnitude of -40 dB can be retained for 
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Figure 6. Estimated nonlinear contributions of the piano soundboard, computed 

according to Eq. ^ and evaluated at three positions of the soundboard. : 

Csaif) for Ai. : Csaif) for A2. : Cs^{f) for A3. An averaging window 

is applied to the amplitude of each spectrum so that each point corresponds to the 
average on a 100 Hz-bandwidth. See text for missing parts. 



the total distortion rate at the ff playing level. 
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3 Experimental modal analysis method 



The experimental study presented in this section aims at estimating the modal 
parameters (modal frequencies, modal dampings and modal shapes) of the 
upright piano soundboard in a wider than usual frequency range (see be- 
low). The terminology and symbols are the followings: when the soundboard 
vibrates freely, each mode corresponds to a time-signal of the generic form 
e~"* cos(27r ft + ip) where / is the modal frequency and a is the modal damp- 
ing, or damping rate. The modal loss-factor is defined by r/ = — - (also: twice 

the damping ratio and the reciprocal of the quality factor). The modal over- 
lap fi is defined as the ratio between the half-power modal bandwidth and the 

, 1 . A/_3dB fv a 
average modal spacmg: /i = — — = — — = — — — . 

A /mode A /mode A / 

The piano was put in a pseudo-anechoic room and excited either mechanically 
and impulsively with an impact hammer or continuously and acoustically with 
a loudspeaker. The former excitation yields the modal shapes but does not 
provide enough energy beyond a certain frequency, depending on the exper- 
imental conditions (here: ^ 500 Hz), hence the acoustical excitation (which 
does not give access to the modal shapes). 

When a loudspeaker was used to excite acoustically the piano, the driving 
signal was adjusted in order to obtain a vibration level similar to the one used 
in the nonlinearity study; the linear contribution of the response was extracted 
prior to modal analysis. The impulsive excitation does not permit to separate 
the linear and nonlinear contributions. Nevertheless, in the light of the results 
presented above and considering the small amplitudes of displacement caused 
by the impacts on the soundboard (typically less than 8 ■ 10~^ m, mostly due 
to a very low-frequency displacement and still less than 1/100 of the board 
thickness) , we considered that the linear approximation was verified up to the 
precision of our measurements. In a way, this was confirmed by the fact that, 
in modal identifications, we never observe a frequency or damping rate that 
was exactly twice or three times that of a lower mode. 

The impact hammer (Kistler - type 9722A) was struck at the nodes of a 
rectangular mesh of 12 x 10 points regularly spaced (Fig. [?]). The motion of 
the soundboard was measured with accelerometers (two B&K 2250A-10 and 
three B&K 4393) at five points in different zones of the board (Fig. 17]). 



For the acoustical excitation of the soundboard, the procedure is the same 
as in Sec. |2] (an exponential swept-sine [50-4000] Hz with a 40 kHz sampling 
frequency, a T = 26 s duration and an assumed N = 4 maximum order of 
nonlinearity). The impulse response of the soundboard is reconstructed by 
the deconvolution technique described in [2] and analysed through a series of 
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Figure 7. Rear view of the upright piano, with the mesh for modal analysis (in red) 
and the locations of the five accelerometers (in black). 



band-pass filters (a typical bank filtering analysis is displayed in Fig. [sjbetween 
550 and 1150 Hz). The cut-off frequencies of the finite-impulse-response (FIR) 
filters were chosen at local minima of the Fourier spectrum of the response. 
If necessary, when there is a doubt on the number of components in one 
frequency-band, two successive overlapping filters were occasionally chosen. 
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Figure 8. Typical bank- filtering analysis of a reconstructed impulse response between 

550 and 1150 Hz (acoustical excitation). )■■■ , : amplitude responses of 

the (slightly overlapping) narrow-band filters. : Fourier spectrum of the impulse 

response at point A2. •: modes estimated by the high-resolution modal analysis 
(modal amplitudes and frequencies). 
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In the piano, the modal overlap factor fi exceeds 30% for all frequencies above 
150 Hz. In other words, almost the whole frequency range of interest is outside 
the low-frequency range where the usual modal analysis technique, based on 
Fourier-analysis, is applicable. With such modal overlap factors, the estimation 
of a damping rates by the width of the peak in the spectrum is not possible 
either. The modal behaviour of the soundboard of the upright piano has been 
investigated by means of a recently published high-resolution modal analysis 
technique [2j which avoids the frequency-resolution limitations of the Fourier 
transform. 



The modal analysis technique used here is well suited for structures made 
of moderately damped materials such as spruce, and for frequencies where 
the modal overlap lies between 30% and 100%. Based on the ESPRIT algo- 
rithm ^Oj, it assumes that the signal is a sum of complex exponentials and 
white noise; it projects the signal onto two subspaces: the subspace spanned by 
the sinusoids (signal subspace) and its supplementary (noise subspace). The 
rotational invariance property of the signal subspace (see |10| for details) is 
used to estimate the modal parameters: frequencies, damping factors and com- 
plex amplitudes. The dimensions of both subspaces must be chosen a priori 
and the quality of the estimation depends significantly on a proper choice for 
these parameters. The best choice for the dimension of the modal subspace 
is the number of complex exponentials actually present in the signal. This 



number (K in Fig. 
modal analysis itse 



9| is twice the number of decaying sinusoids. Prior to the 
f, an estimate of this number is obtained by means of the 
ESTER technique [H] which consists in minimising the error on the rotational 
invariance property of the signal subspace spanned by the sinusoids. The block 
diagram given in Fig. |9] describes the three main steps of the modal analysis 
method: (a) reconstruction of the acceleration impulse response, (b) signal 
conditioning, (c) order detection, (d) determination of modal parameters. 



Below 500 Hz, for each of the 120 x 5 measurements, reconstructing and 
analysing the impulse response as described in Fig. [9] yields the results pre- 



sented in the top frames of Figs. 10 and 11 In order to measure the damp- 
ing with enough precision, it proved necessary to band-filter the impulse re- 
sponses prior to analysis, yielding the results presented in the middle and 



bottom frames of Figs. [TO] and [TT] respectively. Finally, the modal dampings 
were extracted by averaging the results after suppression of the (usually poor) 
estimations in the nodal regions: bottom frame of Fig. 
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Above 550 Hz, the clusters in Fig. [I0[a) are more difficult to identify, owing to 
a too low Signal-to-Noise Ratio (SNR ^ 35 dB): the excitation by an impulse 
hammer is intrinsically limited in frequency since the force exerted by the ham- 
mer is of finite duration (in the order of the time taken by the initial impulse 
to be echoed by the closest discontinuity). Technically, the signal-to-noise ratio 
is not high enough beyond ~ 500 Hz for determining modal parameters with 
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Figure 9. Block diagram of the high-resolution modal analysis method with which 
the modal frequencies fk, modal dampings ak, modal amplitudes ak and phases ipk 
are determined (after [2]). 
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Figure 10. Modal frequencies and dampings in the [0-600] Hz frequency-band, ob- 
tained after an impulse excitation and given by a high-resolution modal analysis, 
(a) Direct analysis, (b) Analysis including a narrow band-pass filtering, (c) Results 
after suppression of the (usually poor) estimations in the nodal regions, (o) : retained 
modal parameters. ( x ) : weighted means of the modal parameters estimated at four 

points of the soundboard with an acoustical excitation (see following section). : 

constant loss-factors r] =1, 2 and 3%. 
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Figure 11. Necessity of the narrow band-pass filtering step. Results for five modes 
in the [230-330] Hz frequency-band before filtering (top diagrams), after filtering 
(bottom diagrams). 

enough precision. In order to extend the estimation of the modal frequencies 
and dampings above 500 Hz, we replaced the impulsive mechanical excitation 
by a continuous acoustical one (Fig. |4]). The impulse responses were extracted 
and the modal frequencies and dampings were estimated by the procedure 
presented in Fig. |9] but the modal shapes could not be determined. 

The results are presented in Sec. [5] and discussed together with the results of 
the numerical modal analysis. 



4 Finite-element model of the piano soundboard 



A finite-element model (FEM) of the soundboard has been written by means of 
the free software CastSM [12j. The soundboard model is that of a rectangular 
plate with two strong bars delimiting two so-called "cut-off corners" in oppo- 
site angles, 11 thin bars (ribs) and the two bridges. Following makers habit, 
the structure is given the shape of a spherical shell, here of radius i? = 43 m, 
corresponding to a value of ^ 8 mm for the so-called crown at the centre of 
the soundboard. This value of the crown is consistent with Conklin's obser- 
vations |13| . corresponds approximately to the plate thickness and represents 
a standard value according to discussions with piano manufacturers. Differ- 
ent manufacturing processes can be employed to realize crowning, all leaving 
residual stresses at the time when the soundboard is fixed in its rim, before 
applying string loading. In the absence of documentation, we have chosen to 
ignore residual stresses here. 
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Figure 12. Mesh of the soundboard: wood panel (black), bars (red), ribs (green) and 
bridges (red). The overall dimensions of the rectangle are 0.91 m x 1.39 m. The 
thickness of the wood panel (without ribs) is 8 mm. The angle between the long-side 
of the soundboard and the "L"-direction of the panel wood (along the grain) is 33°. 

Due to the tension of the strings and to the angle that they form with the plane 
of the soundboard when going through the bridge, a load (called downbearing) 
is exerted on the bridge, in the direction perpendicular to the soundboard. This 
load results in an internal stress and has been included in the numerical model 
as a vertical (transverse) force uniformly distributed along the two bridges. 
This load can be adjusted to some extent by the maker and has been chosen 
here so that the crown at the centre of the soundboard is reduced to one-half 
of its initial (without string loading) value: this is also a standard reduction 
based on piano manufacturers know-how. 

The geometrical data were measured directly on the soundboard (see the cap- 



tion of Fig. 12). All finite elements are triangular thin-shell orthotropic ele- 



ments. The mesh (Fig. 12) has 14267 nodes. The shell thickness is 8 mm; the 
height of the medium-bridge is 3 cm and that of the bass-bridge is 5.5 cm. 
The cut-off bars are 3.3 cm thick. The dimensions of the ribs and the inter-rib 
distances are given in Appendix [C| The boundaries of the soundboard are 
supposed to be clamped and the dynamics is supposed to be conservative. 

All the pieces are made of wood, considered here as an orthotropic material: 
spruce for the rectangular panel and the ribs, fir for the cut-off bars, maple 
for the bridges. For each of these three wood species, four elastic constants are 
necessary to model the chosen finite elements. Pianos are not all made of woods 
with exactly the same characteristics. We have retained the maple character- 
istics given by Haines and the fir characteristics given by Berthaut [15j. 
Although a parametric study falls beyond the scope of this article, we have 
run the FEM simulations with three sets of values for spruce characteristics as 
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given by Haines [H] (Norway spruce) and Berthaut [15] (Sitka spruce). The 
French piano maker Stephen Paulello gave us information on spruce charac- 
teristics that he had observed on ordinary pianos: we refer to it as "mediocre 
spruce". Since the piano that has been chosen here is clearly not a high-end 
model of a well-known brand, the latter wood is plausible; in addition, the 
numerical results also fit best the experimental results, in terms of modal den- 
sity (see Fig. [I6| . The corresponding numerical values and that of p are given 
in Table [T| The "L" (longitudinal) direction refers to that of the grain and 
corresponds to the main axis of orthotropy with the higher elasticity modulus. 
The "R" (radial) direction is across the grain and corresponds to the other 
main axis of orthotropy. The ribs are cut in the "L" direction of their wood. 
They are glued on the panel in the "R" direction of the panel's spruce, corre- 
sponding to Oy. The cut-off bars and the two bridges are also cut in their "L" 
direction. 



In order to obtain the same crown under string loading for each set of spruce 
characteristics, and owing to the possible adjustment of the downbearing by 
the piano makers, the overall static force value was set to 2200 N, 1700 N 
and 1400 N for the Norway spruce, the Sitka spruce and the mediocre spruce 
respectively. 





(GPa) 


(GPa) 


Glr (GPa) 




p (kg m ^) 


Fir 


8.86 


0.54 


1.60 


0.3 


691 


Maple 


10.0 


2.20 


2.0 


0.3 


660 


Norway spruce 


15.80 


0.85 


0.84 


0.3 


440 


Sitka spruce 


11.50 


0.47 


0.5 


0.3 


392 


Mediocre spruce 


8.80 


0.35 


0.5 


0.3 


400 



Table 1 

Density and elastic constants of fir (cut-off bars), maple (bridges) and spruce (ribs 
and panel). The data of the first and fourth lines are given by Berthaut [15|, those 
of the second and third lines by Haines [14J, that of the last line by Paulello. The 
subscripts l and r stand for "longitudinal" and "radial" respectively. The radial 
and longitudinal directions refer to how strips of wood are cut and correspond to 
the "along the grain" and "across the grain" directions respectively. In the geometry 



of the soundboard (Figs. [Tj [7] and 12 ), the x- and y- directions correspond to l and 
R respectively for the spruce panel: = Eji = E,, 



y 



5 Results and discussion 



Results are presented and discussed in terms of modal dampings (Sec. 
first modal shapes, boundary conditions and first modal frequencies (Sec. 



5.1), 



5.2). 
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Finally, the modal density n{f), defined as the number of modes per Hz, or 
reciprocal of the average frequency-interval between two modes is discussed in 
Sec. 15.31 



5.1 Modal dampings 



The modal dampings are reported up to 500 Hz in the bottom frame of Fig. [TO 



and up to 3 kHz in Fig. 13, together with bibliographical results. The observed 



values yield values of the modal overlap fi ranging from around 30% at 150 Hz 
to around 70% at 550 Hz: this explains why the bibliographical results on 
modal dampings, obtained by modal analyses based on the Fourier transform, 
are limited to ^ 500 Hz. The frequency-domain explored here includes mid- 
frequencies, which makes the acoustical excitation technique combined with 
the high-resolution analysis very appealing for modal analyses of musical in- 
struments. 



Except for the first four low- frequency resonances (see Fig. 10 -c), at which the 



energy losses at the rim are probably not negligible compared to those inside 
wood, the modal loss-factors up to around 1200 Hz range from 1% to 3% 
(mean of rj ^ 2.3% for the 55 lowest-frequency estimations). This corresponds 
roughly to what they would be if losses were located in spruce only, where loss 
factors lie commonly between 1 and 3%. 
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Figure 13. Modal frequencies (abscissa) and damping factors (ordinates). (o) : es- 
timations at point A2 of the soundboard. Bibliographical results: o fi6\; (+) |17) : 
(•) [18] . : constant loss-factor curves 77 = 1, 2 and 3%. 



At higher frequencies, modal dampings increase from a mean value of about 
80 s-i below 1200 Hz to about 130 s'^ between 1200 and 1500 Hz. This is 
probably due to a change in the proportion of the energy lost in wood to 
the acoustically radiated energy. Interestingly, Suzuki [TB] noticed on a small 
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grand that " the transition range from less efficient to efficient sound radiation 
is 1-1.6 kHz" . Above 1.8 kHz (or slightly less than that), the loss factors are 
again in the order of the material loss-factors for spruce. 



5.2 Modal shapes, boundary conditions and modal frequencies 



The first (up to 300 Hz) experimental and numerical modal shapes are pre- 



sented in Fig. 14 , together with the modal frequencies. As expected, the numer- 
ical modal frequencies (labeled "Ns", "Ss", "ms" when obtained with Norway 
spruce, Sitka spruce and mediocre spruce respectively) are sensitive to the 
elastic constants and density of wood. Without entering into a discussion of 
parameter or condition sensitivity which is not the focus of this article, one ob- 
serves also that the numerical frequencies are little sensitive to the crown or to 
the string loading condition and that this sensitivity decreases with the modal 
order (Table [2]). These results are quite comparable, in relative magnitude, to 
the numerical results given by Mamou-Mani et al. pi)J: in this reference, see 
Fig. 3 for 6 = 1 (the value that we have chosen for crowning, expressed in the 
units of this referenc^^ and Fig. 4, for a transverse displacement of 4 mmp] 
However, one should consider all these results with care, due to a methodologi- 
cal problem: as explained in Sec. |4| residual stresses induced by manufacturing 
are ignored in our study, as well as in [19j . Since residual stresses are likely to 
induce changes in modal frequencies, comparable to the (small) changes due 
to string loading, the results given in Table [2| as well as the predictions given 
by the model simulated in [19], can only be considered as roughly indicative, 
with respect to crowning and string loading. 



Mode 


(1,1) 


(2,1) 


(3,1) 


(1,2) 


(2,2) 


(4,1) 


(3,2) 


"C" and "L" 


84 


127 


194 


204 


229 


271 


301 


"C", no "L" 


88 


131 


199 


207 


233 


277 


306 


No "C", "L" 


87 


129 


196 


204 


230 


273 


302 


No "C", no "L" 


88 


131 


200 


209 


234 


280 


309 



Table 2 

First numerical modal frequencies, in Hz. Variations with the crown ("C") and 
loading ("L") conditions, with values given in Table [Tj for Norway spruce. 

The modal shapes do not vary significantly for the three kinds of spruce that 
have been simulated (not shown here): this was expected, given the rather 



^ We are skeptical on the practical character of crowning values corresponding to 
6 > 2. 

^ We have reservations on the applicability of string loading that would reduce the 
initial crowning by more than one h unit, here: 8 mm. 
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271 Hz (Ns), 253 Hz (Ss), 234 Hz (ms) 301 Hz (Ns), 279 Hz (Ss), 253 Hz (ms) 
Mode (4,1) Mode (3,2) 

274 Hz 295 Hz 303 Hz 




Figure 14. First modal shapes and modal frequencies: numerical (upper lines) and 
experimental (lower lines). The modal frequencies of the numerical modes are re- 
ported for the Norway spruce (Ns), Sitka spruce (Ss) and mediocre spruce (ms), in 
this succession order. The modal shapes are those given with mediocre spruce. See 
text for the names and the grouping of modes. 



restricted variations of the ratio between Ei^ and E-£i that have been allowed 
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here. The same insensitivity to the loading and crown conditions (hsted in 
Table [21) has been observed. 



The first mode which is observed experimentally, at 81 Hz, has no numerical 
counterpart. Looking at its modal shape, it appears that the boundary condi- 
tions cannot be considered as clamped or hinged all around the soundboard. 
According to visual inspection and therefore, approximately, the soundboard 
seems to be clamped or hinged along the top side, where the case is and 
which is more massive than the rim. Given the nodal line in the middle of the 
soundboard, the other boundary conditions at this particular frequency are 
dominated by inertia rather than elasticity. We therefore labelled this mode 
as "MMCM" (for mass-mass-clamped-mass). Since the boundary conditions 
at / = Hz are necessarily elastic (the piano stays at its place), it follows that 
the first resonance frequency of the boundary mobility, considered here as a 
whole for the sake of simplicity, is below 81 Hz. At least up to a hypothetical 
second resonance frequency, the inertial nature of boundary conditions have 
several consequences. 

(a) The mobility of the boundary decreases with frequency. Therefore, tak- 
ing clamped (or hinged) boundary conditions becomes a better and better 
approximation, as frequency increases. Such constrained boundary conditions 
are generally assumed in the literature, as well as in our FEM simulations. 

(b) Compared to clamped or hinged boundary conditions, the outer nodal 
line moves toward the inside of the soundboard. Indeed, a hint of a nodal line 
can be seen at the bottom left corner of the (1,1) modal shape. However, the 
contrast between the mobility of the boundary and that of the soundboard 
is such that the outer nodal line cannot generally be distinguished from the 
boundary, to the exception of the MMCM and of the (1,1) modes. 



(c) In comparison with clamped or hinged boundary conditions, inertial bound- 
ary conditions raise the modal frequencies (that can also be understood by 
considering the inward shift of the outer nodal line, described above). More- 
over, the relative shift in modal frequency is expected to decrease as frequency 
goes up (asymptotically, the modal frequencies do not depend on the bound- 
ary conditions). This is what is generally observed on the first modes depicted 



in Fig. [T4j the experimental modal frequencies are systematically larger than 
their numerical counterparts. In consequence, choosing the characteristics of 
the soundboard material by fitting the first numerical modal frequencies to 
the experimental ones is not a good idea. The upper modal frequencies must 
be retained instead. This is why we consider that our soundboard is more 
probably cut in "mediocre spruce" than in either Sitka or Norway spruce (see 



Fig. 16). 



(d) The mobility of the boundary is expected to be very low (see above) 
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and to encounter erratic variations along its perimeter. Compared to strictly 
motionless boundary conditions (used in FEM simulations), this is usually a 
cause for doubling some modes (twin modes) as observed for the (2,1), the 
(3,1) and the (4,1) modes which exhibit similar shapes but different modal 
frequencies. Modal families of modes in low frequency have been observed on 
grand pianos, and attributed to the boundary conditions, by Suzuki [16] in 
1986 and Kindel et al. |20j in 1987. Kindel et al. observed up to three or four 
very similar modal shapes, differing mainly in the motion of the edge (the 
rim, also called the case) of the soundboard. Suzuki also noticed that mode 
splitting disappears when several bags of lead shot are put on the rim. He 
named such modes rim resonances, somewhat misleadingly, in our opinion. 

The labelling of the mode (1,2) is somewhat arbitrary. The deformation of the 
numerical mode labelled "(2,2)" is mainly located in the cut-off corner and 
we did not observe any clear experimental correspondent to this mode. One 
reason might be that no accelerometer was put in the cut-off corners during 
the experiment. 




Mode 28 - 722 Hz 



Mode 52 - 1155 Hz 



Mode 79 - 1556 Hz 




Mode 93 - 1728 Hz Mode 119 - 2103 Hz Mode 167 - 2733 Hz 

Figure 15. Examples of modal shapes obtained by the FEM analysis. Below ~ 1 kHz, 
the ribs configuration is not apparent. Above ~ 1 kHz, the vibration is confined 
between the ribs. In the ribbed zone of the soundboard and above that limit, most 
of the modes are localised in one or two areas, extending over a very few inter-rib 
spaces. 



Some of the higher modal shapes of the FEM are represented in Fig. 15 



as computed with the characteristics of mediocre spruce. A transition occurs 
at ~ 1 kHz (see also next section). Above this limit, the ribs confine wave 
propagation, as suggested by Nakamura |21) . observed by Moore et al. [22] 
at 2837 Hz and characterised in terms of wave- numbers by Berthaut [15| : the 
soundboard behaves like a set of more or less coupled structural wave-guides. 
Moreover, localisation seem to occur in this regime of vibration: whereas the 
modal shapes below 1 kHz extend throughout the entire wood panel, one can 
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see that the modal shapes become locahsed above 1 kHz. It must be empha- 
sised that, on plates, periodically spaced stiffeners restrict the possible values 
of the component of the wave-numbers in the direction normal to stiffeners 
(see [23]l24j for example). Localisation of the modes is almost certainly due 
to the non-uniform rib-spacing that is almost always observed in pianos (see 
Tab. C.2 for the particular piano that has been investigated here), as in the 
case of the well-known Anderson localisation of waves in slightly disordered 
structures. 



5.3 Modal density 



Estimations of the modal density of the soundboard are represented in Fig. [16 



up to 3 kHz. They have been obtained as the reciprocal of the moving average 



on six successive estimated modal spacings. In the top frame of Fig. 16, the 
estimation is done independently at four points of measurement (see FigT^for 
the exact locations): the estimated quantity is the apparent modal density at 
that point. The average modal spacing (inverse of the modal density) is around 
22 Hz for the 21 lowest modes, in agreement with comparable low-frequency 
studies (^ 25 Hz for a similar upright piano in ref. [T^ and ~ 22 Hz for a 
baby grand in ref. |16]). 



Although a parametric study is not in the scope of this paper, we present 
also the estimation of the modal density for numerical modes obtained with 
several sets of wood characteristics (see Table [T]), in the bottom frame of 
Fig. [T6j The case of mediocre spruce is also reported in the top frame. The 
frequency evolution n{f) of the modal density reveals two distinct vibratory 
regimes of the structure. 

Below 1.1 kHz, the four experimental sets of results are almost similar. This 
means that the modes which are detected at each measurement point extend 
over the whole soundboard, as confirmed by the experimental and the nu- 



merical modal shapes that are shown in Fig. 14 The modal density increases 
slowly and tends towards a constant value of about 0.06 modes Hz~^: the 
soundboard seems to behave more or less like a homogeneous plate. Consider- 
ing the highest modes in this frequency domain (as explained in the previous 
section), it appears that the simulations with mediocre spruce are those which 
best fit the experimental data. 

The slow rise of n{f) with frequency is characteristic of constrained bound- 
ary conditions [25j . This experimental observation combined with the analysis 
of the lowest modes (see previous section) and with the observation of the 
mounting of the soundboard at its rim lead us to propose the following sim- 
plified scheme for the boundary conditions: the rotational degrees of freedom 
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are blocked and the translational degrees of freedom are massive; as frequency 
increases, this scheme becomes equivalent to clamped boundary conditions. 



For frequencies above 1.1 kHz, n{f), as measured at a given point, decreases 
significantly. Also, n{f) is slightly but consistently different at each of the 
measured locations of the soundboard. In this frequency domain, the FEM 
and the experimental estimations of the modal density differ completely: the 
apparent modal density, estimated at one given point, is roughly the same 
everywhere but not the same as the global modal density given by a numerical 
simulation. This can be explained by the localisation of the vibrations, as 
suggested by Fig. 15 and by the discussion in Sec. 5^ The modal density 
given by the FEM takes into account all the modes of the structure whereas 
the modes detected by one particular accelerometer are only those having a 
significant level where the accelerometer is located. This also explains why the 
modal densities estimated at different locations are different. For example, Ai 
and A5 are near the corners of the soundboard (see Fig. [?]); they belong to 
shorter structural waveguides and "see" less modes than A2 and A3 which 
are located near the centre of the board where the waveguides are longer. 
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10^ 




10^ 

Frequency (Hz) 

Figure 16. Estimations of the modal density in a piano soundboard. Each estimation 
is the reciprocal of the moving average on six successive modal spacings. 
Top frame: experimental determinations (apparent modal density) and numerical 
simulations with mediocre spruce (see below). Modes are measured up to 3 kHz. 
Dots: observed values at points Ai (.), A2 (a), A3 (t), and A5 (*), whose locations 

are given in Fig. [7] H \ \- : numerical modes given by FEM with the wood 

characteristics of mediocre spruce. 

Bottom frame: numerical modes given by FEM with the wood characteristics of 
Norway spruce (□), Sitka spruce (x) and mediocre spruce (H 1 1-)- 

6 Conclusion 



We have applied original techniques to investigate the vibrations of the sound- 
board of an upright piano in playing condition. The nonlinear part of the 
mechanical response to an acoustical excitation could be separated from the 
nonlinear contribution induced by the loudspeaker. At levels of vibration cor- 
responding to ff playing, the nonlinear component of the soundboard vibration 
is ^ 30-50 dB below the linear part, in the [100 - 3500] Hz frequency range. It 
is likely that the main nonlinearity is a consequence of large displacements of 
the soundboard (geometric nonlinearity). If this is true, vibrations at higher 
frequency and corresponding to comparable acoustical levels are not likely to 



25 



generate larger nonlinear components. One may therefore safely retain the or- 
der of magnitude of -40 dB for the distortion rate at the ff nuance, for the 
piano that has been investigated. For larger pianos, such nonlinearities might 
be expected to be even less. This preliminary study shows that a linear model 
is sufficient to predict the main features of the vibro-acoustical behaviour of 
a piano soundboard in playing situations. 

Given the essentially linear response, modal identifications have been per- 
formed between 50 Hz and 3 kHz by means of a novel high-resolution modal 
analysis technique. For the piano, this frequency range belongs mostly to the 
mid-frequency domain since the modal overlap appears to range from 30% at 
150 Hz to 100% at 1 kHz, decreasing down to 60% at 3 kHz. The loss factor ap- 
pears to be maintained between 1 and 3% over several kHz, with a significant 
dispersion but without strong systematic variations. The dispersion might be 
attributed to the different acoustical efficiencies of the different modes of the 
soundboard. Since the loss- factor commonly observed for spruce is about 2%, 
the energy dissipation scheme is likely to be that in which only a small part 
of the power is radiated acoustically at any frequency, thus providing a long 
decay time for each note, at each of its partial frequencies (commonly but in- 
correctly called "harmonics"). On one hand, this raises the hope that using a 
less lossy material than wood may provide a higher sound level together with 
keeping the same decay-rate for each note. However, two other requirements 
for a good tonal quality of the instrument must be kept in mind: spectrum 
regularity (one partial must not behave too differently from the other ones) 
and homogeneity along the tonal range (notes must not differ appreciably, at 
least from their neighbours in pitch). Since the efliciency of the acoustical ra- 
diation of the soundboard is expected to be much more frequency-dependent 
than loss-factors in wood, energy losses that would be caused primarily by 
acoustical radiation must be expected to affect negatively spectrum regularity 
and tonal homogeneity. In other words, the fact that in today's pianos, losses 
seem to be mainly located in wood certainly smoothens frequency-dependency 
of the decay rates, which is favourable for these two timbre qualities. 

The frequency evolution of the estimated modal density of the piano sound- 
board reveals two well-separated vibratory regimes of the structure. Below 
approximately 1 kHz, the modal density and the modal shapes look like those 
of a homogeneous plate. The vibration extends over the whole area of the 
soundboard, including in the so-called "dead-zones". Analysing together the 
modal shapes, the modal frequency and the evolution of the modal density 
in the low-frequency domain suggests that boundary conditions can be con- 
sidered as (a) ruled by inertia for the one or two very low modes and (b) 
constrained in general. The scheme that we propose for boundary conditions 
is that the rotational degrees of freedom are blocked whereas the translational 
degrees of freedom are massive. 
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Above 1 kHz, the soundboard operates as a set of structural wave-guides 
defined by the ribs, as aheady noticed by several authors. The modal shapes 
obtained by FE-modelling confirm this confinement and suggest that modes 
are localised in restricted areas (one or a few inter-rib spaces), due to a slightly 
irregular spacing of the ribs across the soundboard. 

All these observations pave the way for a very synthetic modelling of the 
soundboard vibration up to several kHz. 



Acknowledgements 

This work has been initiated during the PhDs of the first and third authors at 
the LMS, for which they were sponsored by the French Ministry of Research. 
We express our gratitude to Andrei Constantinescu for his precious help on 
the implementation of the FEM of the piano soundboard. We thank Stephen 
Paulello for sharing his knowledge of piano making with us on a number of 
questions. 



A Appendix: Volterra series for modelling weakly nonlinear sys- 
tems 

Volterra series |26] are a means to express the relationship between the input 
e{t) and the output s{t) of any weakly nonlinear system [UfT] as a series of 
multiple convolution integrals: 

s(t)=V/ ... / Vk{Ti,...,Tk)e{t-Ti)...e{t-Tk)dTi...dTk (A.l) 

The functions {vk{Ti, ...,Tk)}k&N* are called Volterra kernels and completely 
characterise the nonlinear system under study. Volterra series are a generali- 
sation of the simple convolution operator used for linear systems. Following the 
idea of linear transfer functions for linear systems, nonlinear transfer functions 
can be obtained by expressing the Volterra kernels in the frequency domain 
via a multidimensional Fourier transform: 

VA;eN* Vk{fi,...,fk)= f Vk{n,...,n)e-'^-^'^\..e-''-f''^'^dn...dn (A.2) 

The association of two weakly nonlinear systems Sp and Sq is now considered, 
as shown in Fig.|5l The families of kernels fk)}keN*, {Gkifi, fk)}keN* 
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and, {Hk{fi, fk)}k£n* fully represent the systems Sp, So and, Sh- In the 
present case, the kernels of Sh can be expressed as functions of the kernels of 
Sp and Sg as follows [26j: 



Hkifl, /fc) — X! X! Pmiifl-, /mi) X • • • 
p=l M'' 

... X i^rn,p(/mi+...+mp_i+l; fk) 

^Gp{fi + . . . + fmii /mi+...+mp_i + l + fn) 



(A.3) 



VA; G N* and with 



m\+...+mp=k 



mi, ... ,mp>l 



For the two first terms, Eq. (A.3) reduces to: 



H2{hj2) = F2(/l,/2)Gl(/l + /2)+Fi(/i)Fi(/2)G2(/l,/2) 



(A.4) 



This proves the intuitive results that the linear transfer function describing 
the linear behaviour of a weakly nonlinear system Sg following another weakly 
nonlinear system iS^ is simply the product of the linear transfer functions of 
those two systems. 



B Appendix: Estimation of the kernels of a cascade of Hammer- 
stein models 



The mathematical foundations of the method used for the estimation of the 
elements of a cascade of Hammerstein models [5|5] are given in this section. 
In such a system, each branch is composed of one nonlinear static polynomial 
element followed by a linear one Pn{t) and the relation between its input e(t) 
and its output s{t) is given by Eq. (B.l), where * denotes the convolution. 



N 

s{t) = Y.(3n*e'\t) (B.l) 

n=l 

To experimentally cover the frequency range on which the system under study 
is to be identified, cosines with time- varying frequencies are used. If e{t) = cos[$(t)] 
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is the input of the cascade of Hammerstein models, the output of the nonhnear 
block e*(t) is rewritten using Chebyshev polynomials as in Eq. (B.2). Details 
of the computation of the Chebyshev matrix C are provided in [8^i9j . 



^ie{l...N} e\t) =cos'[^{t)\ = Y,C{i,k)cos[k^{t)] (B.2) 

fc=0 



When the instantaneous frequency of e{t) is increasing exponentially (from fi 
to /2 in a time interval T), the signal is called "exponential sine sweep". It can 

be shown [8]l9] that by choosing = ^2m7r — — ^ IlJ^ with m G N*, one 
obtains: 
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VA; e N*, cos(A;$(t)) = cos($(t + At/,)) with AtA 



TfYi Iri }\i 
ln(/2//i: 



(B.3) 



which represents another expression of the k term in the linearisation pre- 



sented in Eq. (B.2). 



For any T^-long exponential sine sweep, multiplying the phase by a factor k 



yields the same signal, advanced in time by Atk- Using Eqs. (B.3) and (B.l) 
one obtains: 



TV 



N 



sit) = E 7n * e(t + At„,) with 7„(t) = C{k, n)hk{t) 



(B.4) 



n=l 



k=l 



7„(t) corresponds to the contribution of the different kernels to the n*^ har- 
monic. In order to separately identify each kernel of the cascade of Hammer- 
stein models, a signal e(t), operating as an inverse of e(t) in the convolution 
sense, is needed. The Fourier transform E{f) of the inverse filter e(t) can 
be built in the frequency domain by means of Eq. ( |B.5 ), where E*{f) is the 
complex conjugate of E{f), the Fourier transform of e{t). 



E{f) 



E{f) 



-1 



h/2 -/i]n[/i,/2] 



(/) 



E*{f) 



\E{fW + e{f) 



(B.5) 



e(/) is a frequency-dependent real parameter chosen to be in the bandwidth 
and to have a large value outside of the bandwidth, with a continuous transi- 
tion between the two domains. After convolving the output of the cascade of 



Hammerstein models s{t) given in Eq. (B.4) with y{t), one obtains: 
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N 



(B.6) 



1=1 



Because Atn oc ln{n) and /2 > /i, the higher the order of hnearity n, the 
more advanced is the corresponding 7n(i)- Thus, if Tm is chosen long enough, 
the different 7„(t) do not overlap in time and can be separated by simply 



windowing them in the time domain. Using Eq. (B.7), the family {Pnit)}n€[i,N] 
of the kernels of the cascade of Hammerstein models under study can then be 
fully extracted. 



[Mt)) 



Al 



(B.7) 



(.)^ stands for matrix transposition and Ac is the Chebyshev matrix C defined 
earlier, from which the first column and the first row have been removed. 



C Appendix: Dimensions of the ribs for the FEM 



In order to ease the implementation of the FEM, the geometry of the ribs has 
been simplified: each rib is given a uniform height (or thickness) all along its 



length. In reality, ribs are tapered (Fig. C.l), thus giving less stiffness to the 
soundboard near its edges. In the FEM, the rib thickness is averaged over the 
rib length, thus keeping the same rib mass but not exactly the same rigidity 
since this mechanical property is proportional to h^. The retained geometry 



is given in Tabs. C.l and C.2 




\ 


1 H\ 


1 










(b) 





ft. = 8 I 



Figure C.l. Rib and soundboard geometry, (a) Side view of a rib. (b) Partial front 
view of the ribbed soundboard. 



30 



Rib # 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


Li 


16.0 


18.5 


17.5 


16.5 


16.5 


12.0 


13.0 


15.0 


15.0 


15.0 


10.0 


L2 


8.0 


19.0 


33.5 


48.5 


49.0 


71.5 


59.5 


46.0 


37.0 


21.0 


11.0 


L3 


6.0 


8.0 


11.0 


13.0 


15.0 


15.0 


16.0 


17.0 


16.5 


22.5 


26.5 


a 


2.0 


2.6 


2.8 


2.8 


3.0 


3.2 


3.0 


3.0 


2.8 


2.5 


2.0 


b 


2.3 


2.3 


2.6 


2.6 


2.6 


2.6 


2.6 


2.6 


2.5 


1.8 


1.7 




1.6 


1.8 


2.1 


2.2 


2.2 


2.3 


2.3 


2.2 


2 


1.4 


1.2 



Tabic C.l 

Dimensions of the ribs. Rib ^1 is at the treble end of the soundboard. 6av is the 
averaged thickness of the rib. All values are in cm. 



Rib # 


1-2 


2-3 


3-4 


4-5 


5-6 


6-7 


7-8 


8-9 


9-10 


10-11 


P 


12.7 


12.8 


12.9 


13.5 


13.0 


13.7 


13.0 


13.0 


12.7 


12.8 



Tabic C.2 



Intcr-rib space p (in cm). Note that the spacing docs not sccni to follow a simple 
law and may simply be considered as irregular. 
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